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Abstract 

This manuscript aims to develop and describe gain scheduling control concept for a gas turbine engine which 
drives a variable pitch propeller. An architecture for gain-scheduling control is developed that controls 
the turboshaft engine for large thrust commands in stable fashion with good performance. Fuel flow and 
propeller pitch angle are the two control inputs of the system. New stability proof has been developed for 
gain scheduling control of gas turbine engines using global linearization and LMI techniques. This approach 
guarantees absolute stability of the closed loop gas turbine engines with gain-scheduling controllers. 

Nomenclature 

Ni: Non-dimensional Fan Spool Speed 
N2: Non-dimensional Core Spool Speed 
T: hrust (N) 

TSFC: Thrust Specific Fuel Consumption 
a: Scheduling Parameter 
a: Singular Value 
A: Eigenvalue 



1 Introduction 

The gain-scheduling approach is perhaps one of the most popular nonlinear control design approaches which 
has been widely and successfully applied in fields ranging from aerospace to process control [171 HI] • Gain- 
scheduling, specifically has been used for gas turbine engine control, some of these works are [H [TOl [H 
[2TJ |25j [2j. In general, stability and control of gas turbine engines have been of interest to researchers 
and engineers from a variety of perspectives. Stability of axial fiow fans operating in parallel has been 
investigated in [22]. An application of robust stability analysis tools for uncertain turbine engine systems is 
presented in [Tj. Application of the Linear- Quadratic- Gaussian with Loop-Transfer- Recovery methodology 
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to design of a control system for a simplified turbofan engine model is considered in [8]. A unified robust 
multivariable approach to propulsion control design has been developed in [7 . A simplified scheme for 
scheduling multivariable controllers for robust performance over a wide range of turbofan engine operating 
points is presented in [9]. 

In the previous work by authors fW, controllers developed for single spool and twin spool turboshaft 
system. Those controllers were designed for small transients, and small throttle commands. In this work we 
develop a gain-scheduling control structure for JetCat SPT5 turboshaft engine using the method presented 
in [T8| |T7l [T9j [20] . the controller is designed to be used for entire flight envelope of the twin spool turboshaft 
engine. 

In this manuscript, first a linear representation of the turbofan system dynamics is developed. Then 
control theoretic concepts for gain-scheduling control of this model is presented. The developed controller 
can be used for the entire flight envelope of the engine with guaranteed stability. Finally the simulation 
results for gain scheduling control of a physics-based nonlinear model of the JetCat SPT5 turboshaft engine 
are presented. 

2 Gain Scheduling Control Design 

Consider the nonlinear dynamical system 

y{t)=gP{xnt),um 

where G is the state vector, u G is the control input vector, G is the output vector, /^(.) is 
an n-dimensional, and g^{.) is an /c-dimensional different iable nonlinear vector functions. We want to design 
a feedback control such that y{t) r{t) as t ^ oo, where r{t) G Dr C is the output reference signals 
vector. 

Assume that for each r{t) G Dr^ there is a unique pair (xf,^ie) that depends continuously on r and 
satisfies the equations: 

r = gP{xP,Ue), 

in case of a constant r. is the desired equilibrium point and Ue is the steady-state control that is needed 
to maintain equilibrium at . 

Definition 1. The functions x^{a)^Ue{a), and re{a) define an equilibrium family for the plant ^ on the 
set Q if 

fP{xP{a),Ueia),re{a)) = 0, 
gP{xP{a),Ue{a)) = re(a), a e n. 

Let Q C M^+^ be the region of interest for all possible system state and control vector (x^, u) during the 
system operation, and denote x^* and li*, i G / = 1, 2, as a set of (constant) operating points located at 
some representative (and properly separated) points inside Introduce a set of / regions fti centered at the 
chosen operating points (x^*,ii*), and denote their interiors as l^^o, such that l^jo R^feo = foi" all j ^ /c, 
and |Ji=i = ^- The linear model around each equilibrium point is 

iP = A^(^xP-xr) + Bf{u-u*), 
y = Cf{x-xr)+Df{u-u*)+yt, 
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where the matrices are obtained as fohows 



= 

Cf = 
DP = 



(xf*,«*) 



, V(a;f,w) G 



QfP 

^ I V(a;f,u) G 



(5) 



Here we assume that the common boundary of two regions Qj and belongs to only one of Qj and Qz- 
Note that at each moment, (x^^u) belongs to only one Qi. 

Performing linearizations at a series of trim points gives a linearization family described by 



where 



SxP = AP{a)SxP + BP{a)Su, 
Sy = CP{a)SxP + DP{a)Su. 

SxP = xP — xP{a) 

= y-ye{a), 
Su = 



(6) 



(7) 



— Ue{a)^ Va G ft. 

Gain scheduled controller for plant (|6|, is designed as follows. First, a set of parameter values ai are 
selected, which represent the range of the plant's dynamics, and a linear time-invariant controller for each is 
designed. Then, in between operating points, the controller gains are linearly interpolated such that for all 
frozen values of the parameters, the closed loop system has excellent properties, such as nominal stability 
and robust performance. To guarantee that the closed loop system will retain the feedback properties of the 
frozen-time designs, the scheduling variables should vary slowly withe respect to the system dynamics [18]. 

The parameter a called the scheduling variable in gain scheduling control. Ap (a) , BP (a) ^ CP (a) ^ and 
DP (a) are the parameterized linearization system matrices and (a), Ue{a)^ and ye{<^) are the parameterized 
steady-state system variables, which form the equilibrium manifold of system ([T]). The subscript 'e' stands 
for steady-state throughout this paper. 

Let the controller have the following structure 



= A%a)Sx'' + B%a)[Sy - Sr], 
5u = C%a)Sx'' + D%a)[Sy - 5r], Va G Q. 



(8) 



where 



(9) 



Sx^ = x^ — xl{a) 
Sr = r — re(a), Va G ft. 

A standard realization of the parameterized controller can be written in the following form with the 



reference signal explicitly displayed 













5u 




_ C%a) 


D%a) 


-D%a) 



dx^ 
Sy 
Sr 



Va G n. 



Controller has the general form 



x%t) = nx%t),yit),rit)), 
u{t)=g%x%t),y{t),r{t)), 

with the input and output signals corresponding to the nonhnear plant Q. 
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(10) 



(11) 



The objective in linearization scheduling is that the equilibrium family of the controller (11) match the 
plant equilibrium family, so that the closed loop system maintains suitable trim values, and second the 
linearization family of the controller is the designed family of linear controllers [17 . 

For the equilibrium conditions of plant ([T]) and controller ( pTj ) to match, there must exist a function x^{a) 
such that 

= f%x^^{a),ye{a),re{a)), 



where 

or, 
sr. 

Qr^c I (^c^ (a), (a), re (a)) 5 
-^l(a;c(a),ye(«),r,(a)), Vct S 

So the controller family has the form 

= A''{a)[x'' — x'^{a)] + B'^{a)[y — r], 
u = C''{a)[x'' - xl{a)] + D''{a)[y - r] + We(a), Ma G fi. 



(12) 



A%a) 



(13) 



(14) 



Note that re(Q^) = ye{<^)^ as a result Sy — Sr = y — r. The scheduling parameter a is treated as a parameter 
throughout the design process, and then it becomes a time- varying input signal to the gain-scheduled con- 
troller implementation through the dependence a{t) =p{y{t)). Thus the gain-scheduled controller becomes 



j.c ^ A''{p{y))[x'' - + B''{p{y))[y - r], 

u = C%p{y))[x' - x'My))] + D%p{y))[y - r] + Ue{p{y)). 



(15) 



Linearization of (15) about an equilibrium specified by a gives 

Sx- = A^ia)Sx^ + B^ia)[y - r] - [A-(a)^|M] x [|^(y,(a))(y - r)], 

Su = C^{a)Sx'^ + D%a)[y -r] + - C"=(a)^|^] x [g(y,(a))(y - r) 



(16) 



Comparing this to (10), we see there are additional terms, which we refer to them as hidden coupling 
terms [IT]- 



In order to get rid of these terms we have to design the controller such that 



da 



(17) 



da 



da 



0. 



It is not always easy to come up with solutions to satisfy this condition. In order to make the design 
process easier, we can augment integrators at the plant input, so there is no need for equilibrium control 
value. By augmenting integrators at the plant ([T]) input we obtain 



xP{t) 

ilit) 



p{xP{t)Mt)) 
-vMt) 










_ VcX I _ 



v{t). 



The new controller has the general form 

x'^it) = r{x'^{t),y{t)At)), 
v{t)=g%x%t),y{t),r{t)), 



(18) 



(19) 



with the input and output signals corresponding to the nonhnear plant ([T]). 
Combining (18) and (IT9|) leads to 



' xP{t) ' 




u{t) 








X 



fP{xP{t),u{t)) 

-vMt) 

r{x'={t),gP{xnt),u{t)),r{t)) 








v{t). 



vit) = g'^ix%t),gPixPit),umr{t)), 

gix,r) 



(20) 



Then the closed loop nonlinear system is 



x{t) = f{x{t),r{t)) + Bg{x{t),r{t)), 
= F{x{t),r{t)) 



(21) 



The augmented linear family of systems for (18) becomes 



xP{t) ' 




" AP{a) 


BP{a) 




" SxP ' 


ilit) 







-T]cX I 




Su 





VcX I 



Xaug 

5y=[CP{a),DP{a 




au 



(22) 



SXaug 



Now, the control realization for this system is 

= A'^{a)x'^ + B'^{a)[y — r], 
V = C%a)x'' + D%a)[y - r] Va G 
One of the options for control design is to set controller matrices as follows 

A%a) = A^ = -ej, B%a) = B"" = I, 

which is a kind of PI control, where Ki{a) is the integral gain matrix, and Kp{a) is the proportional 
matrix. 

Hence the control for the augmented system has the final form 







-ej I -I 


V 




Ki{a) Kp{a) -Kp{a 



X 

y 



With these choices for control matrices, the control input is 



X" = / {-Ccx'' + - r)) dr, 

V = -Ki{a)x'' - Kp{a){y - r), Va G n. 

Figure [l} shows schematically how the gain scheduling controller works. 
The linearized closed loop system (22) with controller (23) becomes 



(23) 

(24) 
gain 

(25) 
(26) 



SxP ' 




AP{a) 


BP {a) 







" SxP ' 







Sit 




'ncD%a)CP{a) 


-T]J + D%a)DP{a) 






Su 


+ 




x^ 




B%a)CP{a) 


B%a)DP{a) 






x^ 




-B%a) 


X 










X 







Sr, 



Va G Vt. 



(27) 
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Figure 1: Gain scheduling controller diagram 



For the case where we have simplified output Sy = fe^, (i.e. C^{a) = /, D^{a) = 0) the linearized closed 
loop system (22) with controller (l25|) becomes 



" 5xP ' 






BP (a) 
















5u 




-r]cKp{a) 




-r]cKi{a) 




Su 


+ 


r]cKp{a) 


Sr, Va G ^. 






I 





-ej 




x^ 




-I 




X 






Aai{a) 




X 









(28) 



2.1 Stability Analysis 



In closed loop system (28), let br = 0, and consider the unforced linear time varying system 

i, = Aci{a)x, x{0) = xq, Va G $7, 
Sy = 5xP. 

Assumption 1. The matrix Ad is bounded and Lipschitz continuous as follows 

\\Aci{t)\\<kA, Vi>0, 

\\Ael{t) - AeliT)\\ <LA\\t-T\\, Vi,T>0, 



(29) 



(30) 



Assumption 2. The constant eigenvalues of matrix Aci{y) are uniformly bounded away from the closed 
complex right-half plane for all constant y. 

Theorem 1. Consider system (2^ , under assumptions^ and^ then there exists constants m, X, and e > 
such that if 

\\m\\<^v, VtG[0,T], (31) 

then 

I \x{t) 1 1 < me-^^ I |xo 1 1 , Vt G [0, T] . (32) 

To analyse the stability of the nonlinear closed-loop system, we use a technique known as "global lin- 
earization" developed in [3 . 



Theorem 2. Consider nonlinear system (21\ l, and assume there are a family of equilibrium points (xe,re) 
such that F{xe^re) = 0. Then A^j = ^ e \/x, where S is a polytope, and it is described by a list of its 
vertices, i.e. in the form 

5:= Co{<,...,^SU' (33) 



where A^l,s are obtained by linearizing nonlinear system (21) near equilibrium points (steady state condition), 



and also non- equilibrium points (transient condition). Now, assume their exist a common symmetric positive 
definite matrix P = > such that: 



AX'P<0, Vie {1,2,..., L}. 



(34) 



then system l[21\ ) is absolutely stable. Since by design Aci{a) G \/a, then system l[29\ ) is also stable. 
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If there are no hidden couphng terms involving 8y^ then the design of a stabihzing hnear controher 
family can be assumed to guarantee stability of the linearized closed-loop system in a neighborhood of every 
a ^ Vt. The closed- loop system is not restricted to remain in a neighborhood of any single equilibrium, but 
is assumed to be slowly -varying and to have initial state sufficiently close to some equilibrium in S. Then the 
conclusion is that the closed- loop system remains in a neighborhood of the equilibrium manifold [17 . Using 
results developed in [18^, we can figure out if a system is slowly- varying or not. Here we rewrite theorem 
Theorem 12 from [17]: 



Theorem 3. For plant suppose the gain- scheduled controller (15) is such that there are no hidden 
coupling terms and the eigenvalues of the linearized closed-loop system satisfy Re[X\ < — e < for every 
a ^ Vt. Then given p > there exist positive constants ja and 7 such that the response of the nonlinear 
closed-loop system satisfies the following property. If the exogenous signal \ \r{t)\\ < ji, for t > 0, and if for 
some a e 









- xl{a) - 










Ue{a) 






x%0) 










then 





' xP{t) ' 






u{t) 











xl{p{y{t))) 




<7, 



< p, Vt > 0. 



(35) 



(36) 



2.2 Integration Anti-Windup 

It is desirable to have integral action in the controller since the presence of an integral term eliminates 
steady state error in the controlled variable. Since the allowable values for control inputs are limited, if any 
controller reaches its limit, and error is produced form the difference of the control signal and the actual 
limited signal applied to the plant. This phenomenon is known as integral wind-up. Because of this. Integral 
Wind-Up Protection (IWUP) is used to reduce the effect of the integral term of the controller. 

An approach to IWUP from [131 IS] was adopted for our control problem. The main idea with this approach 
is to decrease the error seen by the integrators. This allows the integrator to increase to an appropriate value 
and decrease the size of the instantaneous change in magnitude when the controller becomes saturated. First, 
the generated control signal to the fuel- metering valve is subtracted from the saturation value. The resulting 
difference is then amplified by an integral feedback gain (IFB) and subtracted from the input to the integrator. 
The IFB is empirically tuned to provide adequate performance. The IFB is not gain scheduled, a constant 
value is sufficient for good performance. 



3 Turboshaft Engine Example 

We apply the developed gain-scheduling controller to a physics-based model of a turboshaft engine driving 
a variable pitch propeller developed in [6l [15] . For a standard day at sea level condition we found five equi- 
librium points for linearizing the dynamics near them. The linearization matrices for these five equilibrium 
points and steady state values of the engine variables and control parameters are: 

• Equilibrium Point 1 (Full Thrust): 

ul = 1.0, = 16 (deg), xl = 1.0, x^ = 0.9524, T* = 255.8685 (A/"), a* = 1.3810, and the matrices 
are 



-5 
3.5 

Kii = 





-2.3 



0.7 0.7 
0.7 0.6 



1.4 
0.63 



, Kpi 





-0.085 

1.2 1.2 
1.2 1.2 



Ci 



(37) 
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• Equilibrium Point 2: 

ul = 0.7, = 16 (deg), xl = 0.8826, = 0.6263, T* = 181.9711 (A/"), a* 
matrices are 



-2.83 
1.20 

Kio -- 



-0.0008 
-2.10 

" 0.6 
0.6 



, B2 



1.14 
0.78 



0.6 
0.5 



1.1 
1.1 





-0.054 

1.1 " 
1.1 



, ^2 = /, 



1.0822, and the 



(38) 



• Equilibrium Point 3 (Cruise): 

ul = 0.4685, = 16 (deg), xl = 0.7264, x^ = 0.5, T* = 70.5125 (TV), a* = 0.8818, and the matrices 



-1.9 0.061 
0.45 -1.1 



Kis = 



0.5 
0.5 



0.5 
0.4 



1.57 
0.3 





-0.023 

1 1 " 
1 1 



(39) 



• Equilibrium Point 4: 

ul = 0.3, u^ = 16 (deg), xl = 0.5327, x^ = 0.3678, T* = 38.155 (N), a* = 0.6473, and the matrices 
are 



-0.85 0.032 
0.32 -0.64 



0.4 0.4 
0.4 0.3 



1.1 
0.17 -0.011 

0.8 0.8 
0.8 0.8 



(40) 



Equilibrium Point 5 (Idle): 



0.145, u^ = 16 (deg), xl = 0.295, x^ = 0.161, T* = 7.317 {N), a* = 0.3361, and the matrices 



are 



A, 



-0.38 -0.0008 
0.26 -0.34 

0.3 0.3 
0.3 0.2 



B,= 



0.7 
0.1 -0.0024 

0.6 0.6 
0.6 0.6 



(41) 



Other controller parameters are set to 

ec = 1, Vc = 3, 



3x7. 



(42) 



To show the stability of the closed loop system, 20 different (10 equilibrium, and 10 non-equilibrium) 
linearizations have been used, to solve inequality (34), in Matlab with the aid of YALMIP [12] and SeDuMi 



[23] packages. The numerical value for the common matrix P is: 



0.639 
0.034 
0.121 
-0.015 
-0.073 
-0.036 



0.035 
0.391 
0.036 
-0.002 
-0.103 
-0.029 



0.121 
0.036 
0.184 
-0.048 
-0.029 
-0.017 



-0.015 
-0.002 
-0.048 
0.130 
0.028 
0.022 



-0.073 
-0.103 
-0.029 
0.028 
0.322 
0.028 



-0.036 
-0.029 
-0.017 
0.022 
0.028 
0.298 
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Figure 2: Norm of closed- loop system matrix 
{\\Aci{t)\\), and its rate of change {\\Aci{t)\\) 



Figure 3: Closed-loop system eigenvalues {X[Aci{a)]) 




Figure 4: Scheduling Parameter {a{t) = \ \x{t)\\) and Figure 5: Norm of measured output of the system 
its rate of change {a{t) = iZr^ii ) (11^(^)11)5 rate of change (||^(t)||) 



To show that the designed gain scheduled controller works properly on JetCat engine we used it to control 
the engine from idle to cruise condition and then again back to idle condition in a stable manner and with 
good performance. Simulation results are shown in figures [2]to|2l] 



Figure jij shows the history of the norm of closed- loop system matrix ||Aci(t)||, and its rate ||Ac^(t)||. As it 
can be seen the figure shows the boundedness of these two variables in accordance with Assumption [T] where 
kA = 7.4539, and La = 1.0106. Figure |3j shows the history of the closed-loop system matrix eigenvalues 
X{Aci}. As it is apparent, all the six eigenvalues remain negative with the time change of the scheduling 
parameter a, and hence satisfies assumption [2] of the stability theorem. 

Figure |4j shows the history of the scheduling parameter which is a = p{y) = \ \y\ \ = ||x||. It is also shows 
the history of the switching function, which is defined based on the norm of the spool speed equilibrium 
values vector. As it is apparent from the plot, engine operated in the vicinity of at least three equilibrium 
points to be able to accelerate from idle to cruise condition. The norm of the scheduling parameter rate 
a{t) = \ ' ^^^^ been plotted. Figure jsj shows the history the norms of the output vector and its 
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rate. This satisfies the condition of Theorem [T] with Cy = 0.0025. Using formulas from [18], we can compute 
m = 496.7476, and A = 0.5271. 
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Figure 6: History of the unforced closed-loop system 
Lyapunov function V{t)^ and its rate of change V{t) 
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Figure 7: Rate of change of reference signals (r) 



Figu re [6} shows the history of the quadratic time varying Lyapunov function of the unforced closed loop 
system ( |29[ ). As it is apparent, V{t) = SX^ P{t)SX ^ is decrescent and bounded from above and below. The 
history of V, shows that it is non-positive for all t > 0, so the exponential stability of the slowly varying 
system (29) with a gains-scheduling controller is guaranteed. Figure [7| shows the rate of change of the 
reference signals for the outputs of the system. The outputs in this simulation are core and spool speed. 
||f|| < 0.15, which corresponds to the assumption of the ||r|| boundedness in the theoremjsj 




time (sec) time (sec) 



Figure 8: Plant states: core and fan spool speeds Figure 9: Controller states 

Figure [Sj shows the history of the plant states which are core and fan spool speeds. Figure [9j shows 
the time histories of the controller states. Figure [Toj shows the phase plot for core and fan spool dynamics. 
Figure pT] shows the time history of the fan and core spool accelerations, i.e. Ni and N2- 

Figures [12] and [13) show the core and fan spool speeds tracking their reference signals. Figure [14] shows 
the history of thrust and it is following its reference command from idle to cruise condition and then back to 
the idle for standard day, sea level condition. Figure [15] shows the control inputs to the augmented system, 
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Figure 10: Core and fan spool speeds vs. core and 
fan spool accelerations 
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Figure 11: Core and fan spools accelerations 
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Figure 12: Output: core spool speed and its reference Figure 13: Output: fan spool speed and its reference 
signal signal 



v{t) = [vi{t)^V2(t)\^ ^ each element corresponding to one of the control inputs to the original system. To 
keep the engine dynamics within the limits of the operation for the available engine model, some limits have 
been defined on the augmented system control input, —0.18 < v(t) < 0.18. This limits will help to keep the 
fuel control input non- negative and also limits the rate of the control, u{t). 

Figure [T6j shows time rate of fuel and prop pitch angle inputs. Figure [T7| shows fuel flow and propeller 
pitch angle histories as control inputs. For better performance and also to keep the engine in the safe range 



of operation limits has been deflned for the augmented control inputs. Figures [18] and [19) show the baseline 
fuel controller integral {Ki{a)) and proportional {Kp{a)) gain matrices histories. These gains have been 
obtained by interpolation using the previously deigned flxed-gain controllers, each one corresponding to a 



equilibrium point of the engine. The numerical values of these gains are mentioned in equations (38) to (40). 

Figure [20) shows JetCat SPT5 turboshaft engine compressor map. In this map the approximate stall line 
and also the operating line for this simulation has been shown. The engine operates in a safe region with 
a big stall margin during its acceleration from idle to cruise and again decelerating back to idle condition. 
Figure 21, shows the histories of turbine temperature, thrust speciflc fuel consumption (TSFC), compressor 
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time (sec) time (sec) 



Figure 16: Rate of change for fuel and prop pitch Figure 17: Fuel and prop pitch angle control inputs 
angle control inputs {u(t)) ('^(^)) 

pressure ratio and corrected air flow rate. 

Using gain-scheduling control technique, this simulation shows the possibility of controlling the linear 
parameter varying model of the turboshaft engine for a large transient throttle. This case study, simulates 
the engine accelerating from idle thrust to the cruise condition and then decelerates back to the idle condition 
in the standard day sea level condition. To tune the controller we need to trade off between the settling 
time and overshoot percentage. For large throttle transients, the existing controller regulates the controlled 
variable, additional limiters, will be developed to protect critical engine variables from exceeding physical 
bounds and to ensure safer operation. 

4 Conclusions 

We developed a gain-scheduling controller with stability guarantees for nonlinear gas turbine engine systems. 
Using global linearization and LMI techniques, we showed guaranteed absolute stability for closed loop gas 
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Figure 18: Controllers integral gain matrix {Ki{a)) Figure 19: Controllers proportional gain matrix 
parameters histories {Kp{a)) parameters histories 
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Figure 20: JetCat SPT5 engine compressor map with Figure 21: Turbine temperature, TSFC, compressor 
operating line overall pressure ratio and air flow rate histories 

turbine engine systems with gain-scheduling controllers. Simulation results presented to show the applica- 
bility of the proposed controller to the nonlinear physics-based JetCat SPT5 turboshaft engine model for 
large transients from idle to cruise condition and vice versa. 
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